ID Potts, Yang-Lee Edges and Chaos. 
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Abstract 

It is known that the (exact) renormalization transformations for the one-dimensional Ising model 
in field can be cast in the form of a logistic map f(x) = Xx(l — x) with A = 4 and x a function of 
the Ising couplings K and h. Remarkably, the line bounding the region of chaotic behaviour in x is 
precisely that defining the Yang-Lee edge singularity in the Ising model. The generalisation of this 
relation between the edge singularity and chaotic behaviour to other models is an open question. 

In this paper we show that the one dimensional g-state Potts model for q > 1 also displays such 
behaviour. A suitable combination of couplings (which reduces to the Ising case for q — 2) can again 
be used to define an x satisfying f(x) = 4a;(l — x). The Yang-Lee zeroes no longer lie on the unit 
circle in the complex z — exp(h) plane for q ^ 2, but their locus is still reproduced by the boundary 
of the chaotic region in the logistic map. 



1 Introduction, Ising 



Yang and Lee [|2[ [3|, later followed by various other authors J4|], provided an important paradigm for 
understanding the nature of phase transitions by looking at the behaviour of spin models in complex 
external fields They observed that the partition function of a system above its critical temperature T c 
was non-zero throughout some neighbourhood of the real axis in the complex external field plane. As 
T —> T c + the endpoints of loci of zeroes moved in to pinch the real axis, signalling the transition. 
When such endpoints occur at non-physical (i.e. complex) external field values they can be considered 
as ordinary critical points with an associated edge critical exponent. This appealing picture was later 
extended by Fisher to temperature driven transitions jgj. 

On any finite graph G n with n vertices the free energy of an Ising-like spin model can be written as 

n 

F(G n , 0, z) = -nh- In J] (z - z k {0)) (1) 
fc=i 

where the fugacity z — exp(h), and h is the (possibly complex) external field. The z k {0) are the Yang-Lee 
zeroes, which in the thermodynamic limit generally condense on curves in the complex z plane. In the 
infinite volume limit n — > oo the free energy per spin is 



F(G oo ,0,z) = -h- d0p(f3,6)Hz-e ie ) (2) 

J —7T 

where p((3, 9) is the density of the zeroes, which can be shown to appear on the unit circle in the complex 
z plane in the Ising case (the Yang-Lee circle theorem). For T > T c or, if one prefers (3 < C , there is a 
gap with p(0, 0) = for \9\ < 9q, and at these edge singularities we have 

P ((3,6)~(o-e y (3) 

which defines the Yang-Lee edge exponent a. This also implies M ~ (6 — Oq)' 7 . Various finite size scaling 
relations relate the Yang-Lee exponent to the other critical exponents [pj and can be used in numerical 
determinations of critical behaviour Q| . 

At first sight there is no apparent reason why the Yang-Lee edge singularity should bear any relation 
to the onset of chaotic behaviour in non-linear maps. The relation between the two in the case of the ID 
Ising model was exposed in Q by considering the renormalization group flow as an example of just such 
a mapping. The partition function for the ID Ising model is given by 



Zjy (K, h) = exp 



N N 
3=1 J = l 



(4) 



where K = and h = -M;, with J the spin coupling and H the external magnetic field, and periodic 
boundary conditions require ctjv+i = o\- The well-known solution to the ID Ising model proceeds by 
expressing Zn(K, h) in terms of the transfer matrix V as Zn = TrV N , where 

Diagonalising V gives the eigenvalues A± = e K jcosh h ± Vsinh 2 /i + e- 4K | and allows us to express 
the partition function as 

Z N = X^ + X N . (6) 

The Yang-Lee zeroes of this partition function in the complex h plane are the N roots of Zn(K, h) = 0, 
which are for real K the solutions of 

Z N = (\ + ) N + (A_y v = & A+ =exp(^)A_ (7) 



1 



where — N < n < N is odd. This gives the N Yang-Lee zeroes h n = i9 r , 



cos 



(H) yW + sinh 2 ^) = l sin(^) cosh(/,„). (8) 



which may also be rewritten as 



5(0 n )= Vl-e-^cos(^) (9) 



In particular we can see that when X — > oo (the zero temperature "transition point" for the ID Ising 
model) the zeroes are uniformly distributed on the unit circle in the complex z plane, as demanded by 
the Yang-Lee theorem. 

So far, so standard. Now note that the recursive renormalization group transformation for the ID 
Ising model can be obtained by demanding that any renormalised couplings K' and h! satisfy 

Zz(K',h')=A N Z N (K,h) (10) 

where A is some renormalization factor. Thinking in terms of a decimation type renormalization scheme 
it is clear that we can satisfy this by taking 

V(K',ti) = A 2 V(K,h) 2 , (11) 

where V is the transfer matrix given in equ.(|^). Viewed geometrically, we are welding two line segments 
together and demanding a suitable rescaling of the couplings, so one copy of the rescalcd transfer matrix 
V(K', h') must serve in place of two copies of the original V(K, h). This leads to the recursion relations 

e 2h' = e 2h cosh(2K + h) 



cosh(21f - h) 



e4K > _ cosh^if ) + cosh(2fr) 

2cosh 2 (/i) ' [ ' 

The crucial observation of fij] was that these recursion relations could be recast by making use of the 
renormalization invariant m — 1 + e 4K sinh 2 (/i) to eliminate h and introducing the variable 

(13) 



(e iK - 1) 



to transform equs.(|l2[) into the logistic map x' = Ax(l — x). This will exhibit chaotic behaviour for 
< x < 1, i.e. if m = 1 + e iK sinh 2 (/i) < which for imaginary external field, h — id, will occur if 
sin 2 (0) > e- iK . 

What has this got to do with Yang-Lee edge singularities? Looking back at equ.([|) we can see that 
the lowest Yang-Lee zero will lie at sin 2 (6>o) = e~ 4K , which is precisely the "boundary of chaos", m < 0, 
in x observed in the renormalization transformation above. One can also identify a gap exponent for 
the chaotic map which is identical to the Yang-Lee exponent a = —1/2 for the ID Ising model pL It is 
natural to ask whether the identification of the onset of chaos in an RG transformation and the Yang-Lee 
edge singularity is a peculiarity of the ID Ising model, or whether other examples of the phenomenon 
exist. In the remainder of the paper we go to answer the question in the affirmative for the ID Potts 
model where one can also obtain the Yang-Lee zeroes explicitly and construct an exact renormalization 
transformation along similar lines to the Ising model. 

2 ID Potts, Yang-Lee 

The partition function for the ID Potts model is given by 



Zn(v,z) =^exp 
{<?} 



N N 



(14) 
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where the S()s are Kronecker deltas and there are now q possible states for each spin a. and we have 
defined y — e K and z — e h for later convenience. We can write down a transfer matrix for this as a q x q 
matrix V(y, z) with q — 2 diagonal elements (y — \)/(yz) 1 ^ q and a 2 x 2 sub-matrix T(y, z) Q [] 



For q = 2 we recover V(K, h) from T(y, z) providing we identify K = 2K, h — 2h. 

The solution proceeds as in the Ising case by writing Zjq{y, z) = trV(y, z) N and diagonalising V 
The dominant eigenvalues A 0i i come from T(y, z) 

Ao.i = \ ((y(l + z) + q - 2) ± y/W Z -*) + 9-2) 2 + (g-l)4z) (yz)-* (16) 
which can be rewritten as 



2 

with 



Ao,i = | ( M- + z± J(z- t\){z t 2 _) ) (yz)-t (17) 



The other q — 2 eigenvalues given by A2 = A3 = . . . = (y — l)(yz)~v play no role in the thermodynamic 
limit. 

The Yang-Lee zeroes z n = e hn , just as for the Ising model, appear as solutions of 

Z N = (\i) N + (Xo) N = o A 1= exp(^)A (19) 
which, upon substituting in the values above for Ao,i, gives 

cos (W) v^" ~ t2+){Zn ~ <2 - } = iSin (2^) (t+i - + Zn) (20) 

which is clearly of the same form as the Ising result in equ.(^|) for general q and reproduces it exactly 
when q = 2 (and we set K = 2K, h — 2h), as it should. The resemblance runs deeper even for general 
q, as noted in M, If we define z = zj (t+tJ) = yz/(y + q — 2) this may be rewritten as 




=-m(^)(l + ^) (21) 

so in the complex z, plane the Yang-Lee zeroes are again uniformly distributed round the unit circle as 
K — > 00 and t + — » 1, t_ 1. 

We now try and pursue the same path with the renormalization group transformation in the Potts 
model as for the Ising model. We once more demand that renormalised couplings y' and z' satisfy 

Zz(y',z')=A N Z N (y,z) (22) 

and again attempt to solve this by taking V(y',z') = A 2 V(y, z) 2 , where V is now the Potts transfer 
matrix. Since only Ao,i are playing any role in the thermodynamic limit we discard the remaining q — 2 
eigenvalues and concentrate our attentions on the sub-matrix T, by demanding T(y', z') — A 2 T(y, z) 2 . It 
is always possible that an infelicitous choice of renormalization transformation could take us outside the 
space of couplings spanned by T, but we shall see that this is not the case here, at least for the symmetric 
choice of transfer matrix that we have made. 



^^We have chosen a slightly different form of the matrix T(y, z) than jSij for convenience in formulating our renormalization 
group transformations. This simply corresponds to different definitions of the ground state energy. 
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We find the following recursion relations 

1 , , „ A 2 



T {y'z') = (yV + z^-l)) 

(y'z')i (yz)i 

l -{y' + q-2) = -^- T ( z ( q -l) + (y + q -2) 2 ) 



(y'z')i (yz)i 

1 -{z') 1 ' 2 = ^l( zy + y + q -2) (23) 



(y'z')i (yz)i 
which can be used to eliminate A giving 

y'z' y 2 z 2 + z(q — 1) 



y> + q - 2 (y + q - 2) 2 + z(q - 1) 

(gQj _ z 1 / 2 (j/0 + j/ + g-2) 
y' + q - 2 ~ (y + g-2) 2 + z( (Z -l)- 
It is then straightforward to show that, as for the Ising model, an invariant exists — in this case 



(24) 



c= fa(l-*) + g-2) 2 (25) 

— so we can use y and C to reduce our recurrence relations to one for y alone. Eliminating z' from ( |2^ ) 
leads to a single recursion relation which can be written as 

c\\ ( ,s z(y{y + q-2) - (q - l)) 2 

v(y +«?-2)-(g-l) = — -2 ■ (26) 

(y (z + l +9-2 



Now we can use (|2^) to write 



C + 4^ + ,-2) ^"" + 1) + < '- 2 '; ,27) 



So we define 

[(C/4)+ g -l] 

[(y-i)d/ + 9-i)] ( j 

and the relation (^6|) is again reduced to the logistic map with the pre-factor 4, 

x'=4x(l-x). (29) 

For C real and positive x is real and negative and so is outside the domain of chaos, but for C < —4(9—1) 
x is positive we have chaos for < x < 1. On the critical line itself, C = — 4(q — 1), we allow ourselves 
the possibility of complex z — \z\e l8 and find from equ.(E5|) that 



(V + Q- 2; _ , .„ 



e w = (M-K" (30) 

y 

where 

cos® = 1-2 / g ,~ 1) 9V (31) 
y{y + q- 2) 

Given our earlier discussion it is no surprise to find that these are precisely the equations defining the 
Yang-Lee edge singularity in the ID Potts model. 

We have thus seen that defining a decimation type renormalization transformation for the ID Potts 
model gives rise to a set of recursion relations which may be reduced using the renormalization invariant 
of equ. (p5|) to a single equation. This may in turn be mapped on to the logistic equation. The boundary 
of the chaotic region for this logistic map is identical to the critical line of the Yang-Lee edge singularity. 
This behaviour is entirely analogous to that seen in the ID Ising model in 

Note that for complex temperatures y, and so x, is complex and defining w = —4(x — |) turns ( |29| ) 
into the Mandelbrot map on the complex plane, 

w' = w 2 - 2. (32) 
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3 Discussion 



The similarity of the Yang-Lee edge singularity for the general q state Potts models in ID and for the 
Ising model was already remarked in The exponent a = —1/2 is identical for all q > 1, and in suitably 
rescaled variables the Yang-Lee zeroes lie on the unit circle at T — for all q > 1. Given this, it is perhaps 
not so surprising that the relation between the renormalization transformations and the Yang-Lee edge 
singularity also survive to general q from the Ising model. Nonetheless, it is intriguing that an exactly 
soluble model has again demonstrated a close relation between the onset of chaos in a renormalization 
map and the Yang-Lee edge singularity. 

We have not discussed the case < q < 1 in this paper, since [|| makes it clear that the correspondence 
with the Ising model is rather less direct for this. For < q < 1 the Yang-Lee edges z + ,z-, lie on the 
positive real z axis for all temperatures. Inside this interval points exist where a third eigenvalue A2 comes 
into play, which would involve the extension of our renormalization transformation to a 3 x 3 matrix. 
For q = 1 on the other hand, the discussion in the previous section suggests that the critical line for the 
logistic map is defined by C — and equ.(|2|) then shows that this translates to z — (y — l)/y and 9 = 
which defines a circle of radius (y — l)/y in the z plane with a "gap angle" 8 = 0. Direct consideration in 
(|] of equs. (|30| , |3l| ) gives the same result for the Yang-Lee zeroes. This demonstrates that the boundary 
of chaos in the renormalization map and the Yang-Lee edge singularity are one and the same for q = 1 
also. 

In a wider context pathologies of approximate real space renormalization transformations within 
phases in higher dimensional models have been discussed in a rigorous manner in Q . Other cases exist 
where recursive non-linear maps are used in the definition of exact partition functions, notably for spin 
models on Bethe lattices (trees). There has been discussion of chaotic effects in such models when 
q < 1 jlO) and the logistic equation has even been observed for a q = 1 state Potts model (related to 
percolation) on a Bethe lattice with co-ordination number 3 plf . Similarly, there have been extensive 
investigations of the effects of frustration in inducing chaotic behaviour in such maps for spin models on 
trees [|l2| [l4|, |l4|, [l6| but no discussion of any Yang-Lee/chaos link in this context. Since results are 
available for the Yang-Lee edge singularity in the Ising model on </> 3 random graphs [jlTj it would be 
worth investigating the correspondence, if any, between chaos and the Yang-Lee edge singularity in this 
case also. This would make it clearer whether the phenomenon reported in |IJ and here was merely a ID 
quirk, or something more general. 

As a final remark we note that the ID Potts models possess a temperature/field (i.e. y, z) duality j|] 

d z+q-1 

y = r~ 

z — 1 

d y+g-i , . 

z = -j=r (33) 

so the relation between the chaotic behaviour discussed here and the endpoints of Yang-Lee (field) zeroes 
can also be couched in terms of Fisher (temperature) zeroes in the dual y D , z D variables. 
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